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BOREAS RSS-7 Regional LAI and FPAR Images From Ten-Day 

AVHRR-LAC Composites 

ling M. Chen, Josef Cihlar 


Summary 

The BOREAS RSS-7 team collected various data sets to develop and validate an algorithm to allow the 
retrieval of the spatial distribution of LAI from remotely sensed images. AVHRR level-4c 10-day 
composite NDVI images produced at CCRS were used to produce images of LAI and the LPAR 
absorbed by plant canopies for the three summer ILCs in 1994 across the BOREAS region. The 
algorithms were developed based on ground measurements and Landsat TM images (Chen and Cihlar, 
1996; Chen, 1996b). The data are stored hi binary image format files. 

Note that some of the data files on the BOREAS CD-ROMs have been compressed using the Gzip 
program. See Section 8.2 for details. 
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1. Data Set Overview 


1.1 Data Set Identification 

BOREAS RSS-07 Regional LAI and FPAR Images From Ten-Day AVHRR-LAC Composites 

1.2 Data Set Introduction 

These Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation (FPAR) images 
were generated in response to the need within the BOReal Ecosystem- Atmosphere Study (BOREAS) 
modeling community for adequate spatial and temporal coverage of LAI and FPAR estimates across 
the BOREAS region. The 10-Day Advanced Very High Resolution Radiometer (AVHRR) composite 
Normalized Difference Vegetation Index (NDVI) images, produced after various processing steps to 
remove artifacts, are very suitable for this purpose. Although algorithms for some cover types for this 
region are not available, these parameter maps were produced as the best estimates at this stage. 
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1.3 Objective/Purpose 

The objective of this project was to provide quantitative spatial distribution of LAI and FPAR for 
the BOREAS region for the 1994 summer Intensive Field Campaigns (IFCs) for the purpose of 
modeling the carbon, water, energy, and trace gas exchange between the boreal ecosystems and the 
atmosphere. 

1.4 Summary of Parameters 

LAI and FPAR absorbed by plant canopies. 

1.5 Discussion 

These LAI and FPAR images are based on AVHRR level-4c NDVI products and a coregistered 
land cover map (Pokrant, 1991). The level-4c product is based on the AVHRR level-4b product, but is 
further processed to remove or mitigate some artifacts caused by the input data or the compositing 
process. The artifacts of concern are atmospheric contamination and bidirectional reflectance effects for 
AVHRR channels 1 and 2, and atmospheric and surface emissivity effects for AVHRR channel 4. The 
processing was carried out at die Canadian Centre for Remote Sensing (CCRS) using software and 
procedures developed in-house (see Section 9 for details). The spatial and temporal coverage of this 
product is identical to that of the level-4b and -4c products. 

1.6 Related Data Sets 

BOREAS RSS-07 LAI, Gap Fraction, and FPAR Data 

BOREAS Level-4c AVHRR-LAC Ten-Day Composite Images: Surface Parameters 
BOREAS Level-4b AVHRR-LAC Ten-Day Composite Images: At-sensor Radiance 


2. Investigator(s) 

2.1 Investigator(s) Name and Title 

Dr. Jing M. Chen 
Dr. Josef Cihlar 

2.2 Title of Investigation 

Retrieval of Boreal Forest Leaf Area hidex From Multiple Scale Remotely Sensed Vegetation 
Indices 

2.3 Contact Information 

Contact 1: 

Jing M. Chen 

Canada Centre for Remote Sensing 

588 Booth Street, 4th Floor 

Ottawa, Ontario 

K1A0Y7 Canada 

(613) 947-1266 

(613) 947-1406 (fax) 

jing.chen@ccrs.nrcan.gc.ca 


Page 2 



Contact 2: 

Josef Cihlar 

Canada Centre for Remote Sensing 

588 Booth Street, 4th Floor 

Ottawa, Ontario 

K1A0Y7 Canada 

(613) 947-1265 

(613) 947-1406 (fax) 

Josef.Cihlar@ccrs.nrcan.gc.ca 

Contact 3: 

Jaime Nickeson 
Raytheon ITSS 
NASA GSFC 
Code 923 

Greenbelt, MD 20771 
(301) 286-3373 

Jaime.Nickeson@gsfc.nasa.gov 


3. Theory of Measurements 

The theories of LAI and FPAR measurements are documented in the Remote Sensing Science 
(RSS)-07 ground LAI and FPAR data document. The theory of AVHRR measurements is given in the 
BOREAS level-4c AVHRR images document, but the relevant portions are copied below. 

The AVHRR is a four- or five-channel scanning radiometer capable of providing global daytime 
and nighttime information about ice, snow, vegetation, clouds, and the sea surface. These data are 
obtained on a daily basis primarily for use in weather analysis and forecasting; however, a variety of 
other applications are possible. The AVHRR-Local Area Coverage (LAC) data collected for the 
BOREAS project were from instruments onboard National Oceanic and Atmospheric Administration 
(NOAA)-9, NOAA-1 1, and NOAA-12 polar orbiting platforms. The radiometers measured emitted and 
reflected radiation in the visible, near-infrared, one middle-infrared, and one or two thermal channels. 

The primary use of each channel and spectral regions and bandwidths on the respective NOAA 
platforms are given in the following tables: 


Wavelength 
Channel (micrometers] 

1* 0.57-0.69 

2 0.72 -0.98 

3 3.52 - 3.95 

4** 10.3-11.4 

5*** 11.4-12.4 


Primary Use 

Cloud and Surface Mapping 

Surface Water Delineation, Vegetation Cover 

Sea Surface Temperature (SST), Nighttime Cloud Mapping 

Surface Temperature, Day/Night Cloud Mapping 

Surface Temperature 


* Channel 1 wavelength for the Television and Infrared Observation Satellite (TIROS)-N flight model 
was 0.5 -0.90 micrometers. 

** For NOAA-7 and-9, channel 4 was 10.3-1 1.3 micrometers. 

*** For TIROS-N, NOAA-6, -8, -10, and- 12 Channel 5 duplicates Channel 4. 
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The wavelength ranges at 50% relative spectral response (in micrometers) of the bands for the 
platform- specific instruments are: 


Band 

NOAA- 9 

NOAA- 

■11 

NOAA- 

■12 

NOAA- 

■14 

1 

0 . 570 - 0 . 699 

0.572 - 

0 . 698 

0.571 - 

0 . 684 

0 . 570 - 

0.699 

2 

0 .714 - 0 . 983 

0.716 - 

0 . 985 

0.724 - 

0 . 984 

0 .714 - 

0 . 983 

3 

3 . 525 - 3 . 931 

3.536 - 

3 . 935 

3.554 - 

3 . 950 

3.525 - 

3 . 931 

4 

10.334 - 11.252 

10.338 - 

11.287 

10.601 - 

11 . 445 

10.330 - 

11.250 

5 

11.395 - 12.342 

11.408 - 

12.386 

10.601 - 

11 . 445 

11.390 - 

12.340 


The AVHRR is capable of operating in both real-time and recorded modes. Direct readout data were 
transmitted to ground stations of the automatic picture transmission (APT) class at low resolution (4 x 
4 km) and to ground stations of the high-resolution picture transmission (HRPT) class at high 
resolution (lxl km). AVHRR HRPT data were received for the BOREAS region by the OCRS 
Prince Albert Satellite Station (PASS). The LAI and FPAR images in 1994 in this data set were 
produced using the NOAA-1 1 sensor. 


4. Equipment 


4.1 Sensor/Instrument Description 

The AVHRR is a cross-track scanning system featuring one visible, one near-infrared, one 
middle-infrared, and two thennal channels. The analog data output from the sensors are digitized 
onboard the satellite at a rate of 39,936 samples per second per channel. Each sample step corresponds 
to an angle of seamier rotation of 0.95 milliradians. At this sampling rate, there are 1 .362 samples per 
instantaneous field of view (IFOV). A total of 2,048 samples is obtained per channel per Earth scan, 
which spans an angle of +/-55.4 degrees from nadir. 

4.1.1 Collection Environment 

The NOAA satellites orbit Earth at an altitude of 833 km. From this space platform, the data are 
transmitted to a ground receiving station. 

4.1.2 Source/Platform 

Launch and available dates for the TIROS-N series of satellites from CCRS are: 


Satellite 

Launch Date 

Date 1 

Range 

TIROS-N 

13-Oct-l 978 

1 9-Oct-l 978 

to 

30-Jan-1980 

NOAA- 6 

27- Jun-1 97 9 

21-Aug-1984 

to 

23-Jan-1986 

NOAA-B 

29-May-1980 

Failed to achieve orbit 

NOAA- 7 

23-Jun-l 981 

24- Jul-1 983 

to 

30-Dec-l 98 4 

NOAA- 8 

2 8-Mar-l 983 

24-Jul-1983 

to 

1 3 -Aug-1 98 5 

NOAA- 9 

12-Dec-l 984 

1 6-Sep-l 98 5 

to 

1 9-Mar-l 9 95 

NOAA- 1 0 

17-Sep-l 986 

ll-Oct-1986 

to 

1 5-Nov-l 9 93 

NOAA- 1 1 

24-Sep-1988 

28-Jun-1989 

to 

1 3-Sep-l 9 94 

NOAA- 12 

14-May-l 991 

1 1-Aug-l 9 93 

to 

present 

NOAA- 1 4 

30-Dec-l 9 94 

15-May-l 9 95 

to 

present 


4.1.3 Source/Platform Mission Objectives 

The AVHRR is designed for multispectral analysis of meteorologic, oceanographic, and 
hydrologic parameters. The objective of the instrument is to provide radiance data for investigation of 
clouds, land water boundaries, snow and ice extent, ice or snow melt inception, day and night cloud 
distribution, temperatures of radiating surfaces, and SST. It is an integral member of the payload on 
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the advanced TIROS-N spacecraft and its successors in the NOAA series, and as such contributes data 
required to meet a number of operational and research-oriented meteorological objectives. 

4.1.4 Key Variables 

Emitted radiation and reflected radiation. 

4.1.5 Principles of Operation 

The AVHRR is a four- or five-channel scanning radiometer that detects emitted and reflected 
radiation from Earth in the visible, near-infrared, middle-infrared, and thermal-infrared regions of the 
electromagnetic spectrum. A fifth channel was added to the follow-on instrument designated 
AVHRR/2 and flown on NOAA-7, -9, -1 1, and -14 to improve the correction for atmospheric water 
vapor. 

4.1.6 Sensor/Instrument Measurement Geometry 

The AVHRR is a cross-track scanning system. The IFOV of each sensor is approximately 1 .4 
milliradians, giving a spatial resolution of 1.1 km at the satellite subpoint. There is about a 36 percent 
overlap between IFOVs (1.362 samples per IFOV). The scanning rate of the AVHRR is six scans per 
second, and each scan spans an angle of +/-55.4 degrees from the nadir. 

4.1.7 Manufacturer of Sensor/Instrument 

ITT Aerospace 

4.2 Calibration 

No in-flight visible channel calibration is performed. 

4.2.1 Specifications 

IFOV 

RESOLUTION 
ALTITUDE 
SCAN RATE 
SCAN RANGE 
SAMPLES /SCAN 

4. 2. 1.1 Tolerance 

A signal-to-noise ratio of 3: 1 at 0.5-percent albedo. 

4.2.2 Frequency of Calibration 

The Naval Research Faboratory's (NRF's) TIROS-N calibration overlay performs the calibration 
on blocks of telemetry data. For FAC/HRPT acquisitions, a block consists of 20 scan lines. 
Calibration begins by reading the calibration parameters into memory. For each scan line of telemetry 
in a block, the following process takes place: 

• Telemetry data are extracted and unpacked. 

• Ramp calibration data for each of the five channels are decommutated. 

4.2.3 Other Calibration Information 

None given. 


1 . 4 mrad 
1 . 1 km 
833 km 

360 scans/min (1.362 samples per IFOV) 
-55.4 to 55.4 degrees 

2,048 samples per channel per Earth scan 
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5. Data Acquisition Methods 


The BOREAS level-4c AVHRR-LAC images were acquired through the OCRS. Some radiometric and 
geometric corrections are applied to produce the imagery in a spatially corrected form (Lambert 
Conformal Conic (LCC) projection). A full level-4c AVHRR-LAC image contains approximately 
1,200 pixels hi each of approximately 1,200 lines. Before any geometric corrections, the ground 
resolution ranges from 1.1 km at nadir to 2.5 km x 6.8 km at the scanning extremes. Each pixel value 
is stored in a 2-byte field starting with level-4b products. The level-4c images were processed through 
software developed at CCRS. The raw data are available from the CCRS PASS. 

6. Observations 


6.1 Data Notes 

None. 

6.2 Field Notes 

None. 


7. Data Description 


7.1 Spatial Characteristics 


7.1.1 Spatial Coverage 

The spatial characteristics for these LAI and FPAR images are the same as for the BOREAS 
level-4c AVHRR-LAC imagery, which covers the entire 1,000-km by 1,000-km BOREAS region. 
This contains the Northern Study Area (NSA), the Southern Study Area (SSA), and the transect region 
between the SSA and NSA. 

The North American Datum of 1983 (NAD83) coordinates for the comers of the LAI and FPAR 
images are: 

Latitude Longitude 


Northwest (1,1) 
Northeast (1,1200) 
Southwest (1200,1) 
Southeast (1200,1200) 


59 . 36395°N 
61 . 01294°N 
48 . 83387 °N 
50 . 02993°N 


115 ,40859°W 
93 . 28553°W 
110 ,25229°W 
93 . 7 3857 °W 


The northwest comer has a distance ( 1 109.76 km west, 7900.04 km north) from the origin (95°W and 
0°N) of the LCC coordinate. The pixel size is exactly 1 km. 

The NAD83 corner coordinates of the BOREAS region are: 


Northwest 

Northeast 

Southwest 

Southeast 


Latitude 


59 . 97 9°N 
58 . 84 4 °N 
51 . 000°N 
50 . 0 8 9 °N 


Longitude 


111 . 000°W 
93 . 502°W 
111 . 000°W 
96 . 97 0 °W 
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The NAD83 corner coordinates of the SSA are: 


Northwest 

Northeast 

Southwest 

Southeast 


Latitude 


54 . 3 1 9 °N 
54 . 223°N 
53 . 51 3 °N 
53 . 4 1 9 °N 


Longitude 


106 .227°W 
104 . 236°W 
106 . 320 °W 
104 . 368°W 


The NAD83 corner coordinates of the NSA are: 


Northwest 

Northeast 

Southwest 

Southeast 


Latitude 


56 . 2 4 9 °N 
56 . 083°N 
55 . 542°N 
55 . 37 9°N 


Longitude 


98 . 824 °W 
97 . 241 °W 
99 . 0 45 °W 
97 . 489°W 


7.1.2 Spatial Coverage Map 

Not available. 

7.1.3 Spatial Resolution 

Before any geometric corrections, the spatial resolution varies from 1.1 km at nadir to 
approximately 2.5 x 6.8 km at the extreme edges of the scan. The level-4b composite AVHRR-LAC 
images have had geometric corrections applied so that the pixel size is 1 km in all bands. Only pixels 
with view zenith angles 57 degrees or less are used in level-4c product. The LAI and FPAR images 
have the same 1 .0-km spatial resolution. 

7.1.4 Projection 

The coordinate system is the LCC with the two standard parallels at 49 degrees north and 77 
degrees north, respectively, and the meridian at 95 degrees west. 

7.1.5 Grid Description 

The LAI and FPAR images are the same as the level-4 images, which are projected into the LCC 
projection described in Section 7. 1.4 at a resolution of 1 .0 km per pixel (grid cell) in both the X and Y 
directions. 

7.2 Temporal Characteristics 

7.2.1 Temporal Coverage 

The overall time period of AVHRR data acquisition in 1994 was from 09-Apr through 10-Sep. 
Ten-day composite images of NDVI for this period were produced at level-4c. Three NDVI composite 
images corresponding to the IFC periods were used for LAI and FPAR calculation (see Section 7.2.2). 

7.2.2 Temporal Coverage Map 

The 1994 compositing periods from the level-4c data set used for the LAI and FPAR parameter 
images were the following 10-day periods, one during each IFC: 

May 21-31 
July 21-31 
September 1-10 
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7.2.3 Temporal Resolution 

The daily images are composited into nominally cloud-free images over 10-day periods. One 
10-day compositing period was selected to represent each 1994 IFC to create the three LAI and FPAR 
images for this data set. 

7.3 Data Characteristics 

7.3.1 Parameter/Variable 

LAI 

FPAR 

7.3.2 Variable Description/Definition 

LAI is defined as one half the total leaf area per unit ground surface area. For four-sided spruce 
needles, two sides are included. To derive LAI from the digital numbers (DN) in the image, use the 
following formula: 

LAI = (DN-1)/10 

Green FPAR is the fraction of incident PAR that is absorbed by the green leaves in the canopy. It 
excludes the fraction reflected back to space and the fraction absorbed by the background (moss, soil 
and understory in forest, and soil in the crops), but it includes the small fraction that is reflected by the 
background and absorbed by the green leaves on the way back to space. To derive FPAR from the 
DNs in the image, use the following formula: 

FPAR = (DN-1)/100 

7.3.3 Unit of Measurement 

Both LAI and FPAR are unitless; for LAI, it can be expressed as m 2 of leaf area/m 2 of ground 
surface area. 

7.3.4 Data Source 

The level-4c AVF1RR data used to create this product was processed and provided by the CCRS. 

7.3.5 Data Range 

Both the LAI and FPAR images are 8-bit images; i.e., the DNs vary from 0 to 255. 

7.4 Sample Data Record 

Not applicable to image data. 


8. Data Organization 


8.1 Data Granularity 

The smallest unit of data tracked by the BOREAS Information System (BORIS) is the entire set of 
AVHRR parameter image files, the readme fde, and the coordinate conversion code. 

8.2 Data Format(s) 
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8.2.1 Uncompressed Data Files 

This data product consists of the following data files: 


logical 

file description record size 


1 Header /Readme 80 

2 FPAR_AVHRR_IFC1_94 . IMG 1200 

3 FPAR_AVHRR_IFC2_94 . IMG 1200 

4 FPAR_AVHRR_IFC3_94 . IMG 1200 

5 LAI_AVHRR_IFC1_94 . IMG 1200 

6 LAI_AVHRR_IFC2_94 . IMG 1200 

7 LAI_AVHRR_IFC3_94 . IMG 1200 

8 lcc 1200. f 80 


The header file contains descriptive information about the files in the product. 

Each of the composite images contain 1,200 8-bit ( 1-byte) pixels in each of the 1,200 lines. The 
image file fonnat is the same as that of any band sequential (BSQ) image file. 

The FORTRAN program, "lcc_1200.f", is provided for the calculation of longitude and latitude for 
any given pixel and line in the image. When using the software, you need to input the pixel and line 
numbers relative to the top-left pixel, which is ( 1,1). A version of the lcc_1200 code written in C is 
available on request from Jing M. Chen at CCRS. 

8.2.2 Compressed CD-ROM Files 

On the BOREAS CD-ROMs, the ASCII header file for this image is stored as ASCII text; 
however, the image files have been compressed with the Gzip (GNU zip) compression program 
(file_name,gz). These data have been compressed using gzip version 1.2.4 and the high compression 
(-9) option (Copyright (C) 1992-1993 Jean-loup Gailly). Gzip uses the Lempel-Ziv algorithm (Welch, 
1994) also used in the zip and PKZIP programs. The compressed files may be uncompressed using 
gzip (with the -d option) or gunzip. Gzip is available from many Web sites (for example, the ftp site 
prep.ai.mit.edu/pub/gnu/gzip-*.*) for a variety of operating systems in both executable and source 
code fonn. Versions of the decompression software for various systems are included on the 
CD-ROMs. 


9. Data Manipulations 


9.1 Formulae 

9.1.1 Derivation Techniques and Algorithms 

9. 1.1.1 LAI Maps 

The NDVI was first converted to the simple ratio (SR) using 
SR = (1+NDVI)/(1-NDVI) 

hi the calculation of LAI from SR, different algorithms were used for the following 10 land cover 
types: water, mixed wood, deciduous, conifer, transitional forest, tundra, barren lands, cropland, 
rangeland/pasture, and built-up areas. For water, barren lands, and built-up areas, LAI is set to zero. 
An algorithm for the cropland is formulated using available NDVI-LAI relationships in the literature 
(Holben et ah, 1980; Gardner and Blad, 1986; Aase et ah, 1986; Wiegand et al., 1992; Li et al., 
1993). The algorithm for boreal conifer forests is derived from our recent work (Chen and Cihlar, 
1996). The methodology for ground truth LAI measurements is described in Chen and Cihlar (1995a 
and 1995b) and Chen (1996a). The algorithm for deciduous forests is also based on our own work, 
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but is less accurate because of the vegetation dynamics and insufficient field data for the seasonal 
coverage. Mixed wood and transitional forests are considered as the intermediate case between conifer 
and deciduous forests. Tundra and rangeland/pasture are treated the same as cropland because of lack 
of data. Summer images are less useful for determining the overstory NDVI in conifer canopies 
because of the understory effect (Chen and Cihlar, 1996); therefore, the IFC-2 and IFC-3 LAI values 
for conifer species were taken to be 1. 12 and 1.05 times the IFC-1 values. A ceiling of LAI values was 
used for each IFC to remove 'speckles' (unreasonably large LAI values caused by noise in NDVI 
data). These ceilings were 5.5, 6.0, and 5.7 for IFC-1, IFC-2, and IFC-3, respectively. 

The algorithms for boreal forests were developed by correlating SR from Landsat-5 Thematic 
Mapper (TM) with ground-based LAI measurements. In the application of the algorithm to AVHRR 
images, a factor of 1. 10 was used to increase the NDVI of AVHRR to account for sensor differences. 

Caution should be taken when using these LAI values for the calculation of FPAR because some 
models only use LAI inputs, while FPAR is calculated from LAI. Because of the canopy architectural 
difference between the species, it is suggested that the LAI values should be reduced by 50% for 
conifers, 25% for deciduous and mixed wood, and 10% for agricultural crops and tundra in the 
calculation of radiation interception using Beer's Law. 

The equations used in the LAI algorithm are: 


IFC-1 and IFC-3 : 

LAI=0 .594 (SR- 2 .781) 

For 

LAI=0 . 475 

(SR- 2 .781) 

For 

LAI=0 .792 

(SR- 2 .781) 

For 

LAI=1 .188 

(SR- 2 .781) 

For 

LAI=0 . 325 

(SR-1 .5) 

For 

LAI=0 . 325 

(SR-1 .5) 

For 

LAI=0 . 325 

(SR-1 .5) 

For 

LAI = 0 


For 

IFC-2 : 

LAI=0 .493 

(SR- 3 . 637) 

For 

LAI=0 .394 

(SR- 3 . 637) 

For 

LAI=0 . 657 

(SR- 3 . 637) 

For 


LAI=1.12*LAI in IFC-1 For 
LAI=0 . 325 (SR-1 . 5) For 
LAI=0 . 325 (SR-1 . 5) For 
LAI=0 . 325 (SR-1 . 5) For 
LAI=0 For 


mixed wood 
deciduous forest 
transitional forest 
conifer forest 
tundra 
cropland 

rangeland and pasture 
water, barren, and built-up 

mixed wood 
deciduous forest 
transitional forest 
conifer forest 
tundra 
cropland 

rangeland and pasture 
water, barren, and built-up 


areas 


areas 


As shown above, the background SR values are all the same because of the lack of field data. The 
background values were obtained through regression for the conifer species. The IFC-2 LAI values 
were calculated from IFC-1 values to minimize the background effect (Chen and Cihlar, 1996). 

9. 1.1.2 Daily Green FPAR Maps 

Green FPAR refers to the fraction absorbed by green leaves only after the removal of the 
contribution of the supporting woody material to the PAR absorption. The instantaneous green FPAR 
is integrated over the day with a weight equal to the cosine of the solar zenith angle to obtain the daily 
green FPAR presented in the map. The daily green FPAR can be used as a parameter to convert the 
daily absorbed PAR to daily total incident PAR. The dates for IFC-1, IFC-2, and IFC-3 are the same 
as those for LAI maps. 

Similar to the calculation of LAI, different algorithms for the daily green FPAR were used for the 
10 land cover types. FPAR for water, barren lands, and built-up areas is assumed to be zero. An 
algorithm for the cropland is formulated using measurements from crop fields (Chen, unpublished). 
The algorithm is similar to that of Asrar et al. (1984). The algorithms for boreal conifer forests are 
published in Chen (1996b). For deciduous cover, the algorithm is formulated using eight sites in the 
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midsummer (Chen, unpublished). The algorithms are linear relationships between SR and FPAR. 
Tundra and rangeland/pasture are treated the same as cropland because of lack of data. 

The values for boreal forests within the map are the most reliable. FPAR of the overstory is 
calculated from the downwelling and upwelling PAR measurements at two levels: above and below the 
canopy, for over 30 stands. The below-canopy measurements were made using the Tracing Radiation 
and Architecture of Canopies (TRAC) instrument (BOREAS RSS-07 Ground Measurements of LAI 
and FPAR). Because the downwelling PAR through the overstory and upwelling PAR reflected from 
the forest floor are highly variable, average values were obtained from closely spaced measurements 
over long transects (50-340 m). Simultaneous measurements of the downwelling and upwelling PAR 
above the forests were made on micrometeorological towers established for BOREAS. The data were 
collected at the beginning, middle, and end of the growing season to consider the effect of seasonal 
vegetation dynamics and the change in solar zenith angle (Chen, 1996b). 

The equations used in the FPAR algorithm are: 


I FC - 1 and I FC - 3 : 

FPAR=0 . 170 (SR- 2 . 044) 
FPAR=0 . 147 (SR- 2 . 044) 
FPAR=0 . 176 (SR- 2 . 044) 
FPAR=0 . 221 (SR- 2 . 044) 
FPAR=0 .138 (SR-1 .5) 
FPAR=0 .138 (SR-1 .5) 
FPAR=0 .138 (SR-1 .5) 
FPAR=0 


For mixed wood 
For deciduous forest 
For transitional forest 
For conifer forest 
For tundra 
For cropland 

For rangeland and pasture 

For water, barren, and built-up areas 


I FC - 1 and I FC - 3 : 

FPAR=0 . 147 (SR- 3 . 074) 
FPAR=0 . 127 (SR- 3 . 074) 
FPAR=0 . 154 (SR- 3 . 074) 
FPAR=1 . 05*FPAR in IFC-1 
FPAR=0 .138 (SR-1 . 5) 
FPAR=0 .138 (SR-1 . 5) 
FPAR=0 .138 (SR-1 . 5) 
FPAR=0 


For mixed wood 
For deciduous forest 
For transitional forest 
For conifer forest 
For tundra 
For cropland 

For rangeland and pasture 

For water, barren, and built-up areas 


In the equations the background SR changes between IFCs. The background values are different 
from the LAI equations because the SR-FPAR relationship is not strictly linear. These values were 
found from the regressions for the conifer species. A refinement of the algorithm would be to use the 
Modified SR (MSR, Chen 1996c), which is more linearly related to FPAR. IFC-2 LAI values for 
conifers are directly related to IFC-1 values because the background effect in IFC-1 is smaller (Chen 
and Cihlar, 1996). This assumes the normal seasonal growth situation. 

Because of lack of complete field data for all cover types, the most reliable data are for conifer, 
transitional forest, and cropland. The images with quality quote 04 (on a scale of 05) will be made 
available upon request. Please notify the contacts listed in Section 2.3 regarding any problems you 
encounter. 


9.2 Data Processing Sequence 

The processing steps from NDVI to LAI and FPAR are summarized above (Section 9.1); the 
processing steps leading to NDVI are described in the document for BOREAS level-4c AVHRR-LAC 
images. The relevant descriptions are copied below for convenience. 

The level-4c processing sequence is called Atmosphere, Bidirectional and Contamination 
Corrections of CCRS (ABC3) and is described in more detail by Cihlar et al. (1997a, 1997b). 
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9.2.1 Processing Steps 


Step 1: Top-of-the-Atmosphere (TOA) reflectance 

TOA reflectance for channel 1 or 2 is calculated from the corrected TOA radiance, L*(new), with 
the formula given by Teillet (1992). Values of gain G and offset O were calculated with consideration 
of postlaunch sensor degradation (Teillet and Holben, 1994). 

Step 2: Atmospheric correction of AVHRR channels 1 and 2 

The algorithm Simplified Method for Atmospheric Correction (SMAC) (Rahman and Dedieu, 

1994) was used in the processing. The processing was carried out assuming a water content of 2.3 
g/cm 2 and ozone content of 0.319 cm-atm. A constant value of 0.05 was used for optical depth at 550 
nm. The corrections were computed on a pixel basis using solar zenith, view zenith, and relative 
azimuth channels. 

Step 3: Identification of contaminated pixels 

A new procedure was developed to identify the 'contaminated' pixels; i.e., pixels where the 
surface vegetation or soil signal is obscured (Cililar, 1996). The procedure, dubbed Cloud Elimination 
from Composites using Albedo and NDVI Trend (CECANT), is based on the high sensitivity of NDVI 
to the presence of clouds, aerosol, and snow. Three features of the annual surface reflectance trend are 
used: the high contrast between the albedo (represented by AVHRR channel 1) of land, especially 
when fully covered by green vegetation, and clouds or snow/ice; the average NDVI value (expected 
value for that pixel and compositing period); and the monotonic trend in NDVI. Four thresholds are 
required in CECANT to identify a partially contaminated pixel (i,j,t) where i and j are pixel coordinates 
and t is the compositing period: 


Cl(t): 

The maximum channel 1 reflectance of a clear- sky, 
snow- or ice-free land pixel hi the data set. 

Rmin(t): 

The maximum acceptable deviation of the measured 
value NDVI(i,j,t) below the estimated NDVIa(i,j,t). 

Rmax(t): 

The maximum acceptable deviation of the measured 
value NDVI(i,j,t) above the estimated NDVIa(i,j,t). 

Zmax(t): 

The maximum acceptable deviation of the measured 
value NDVI(i,j,t) above the estimated NDVImax(i,j,t). 


NDVImax(i,j,t) and NDVIa(i,j,t) were calculated using the Fourier- Adjustment, Solar Zenith 
Angle Corrected, Interpolated, Reconstructed (FASIR) model of Sellers et al. (1994), which 
approximates the seasonal NDVI curve with a third-order Fourier transform. Before the computation, 
missing NDVI values between first and last measurements were replaced through linear interpolation 
after finding the seasonal peak for each pixel, using the rationale and algorithm of Cililar and Howarth 
(1994). A constant value of 0.30 was used for Cl. The upper and lower limits for R and Z were 
determined separately for each composite period using R and Z histograms (Cililar, 1996). Using these 
thresholds, a cloud mask was prepared for each composite period. 

Step 4: Corrections for bidirectional reflectance effects in channels 1 and 2 

The model of Roujean et al. (1992) as modified by Wu et al. (1995) was used to characterize the 
seasonal bidirectional reflectance function for each cover type, Fand cover-dependent model 
coefficients were derived (Fi et al., 1995) using a map of Canada with pixel size of 1 km prepared with 
AVHRR data (Pokrant, 1991). Only cloud-free pixels were included hi the derivation of the model 
coefficients, and no bidirectional corrections for snow- or ice-covered areas were made. The resulting 
models were used to compute channel 1 and 2 reflectance for view zenith of 0 degrees and solar zenith 
of 45 degrees. 
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Step 5: Replacement of contaminated pixels for AVHRR channels 1 and 2 

Two cases were recognized: pixels contaminated either during or at the end of the growing season. 
For pixels contaminated during the growing season, the new values were found through linear 
interpolation for both channels 1 and 2. At the end of the growing season, it was assumed that the 
annual trajectory for individual channels as well as for NDVI could be approximated by a 
second-degree polynomial. The polynomial was fitted to the plot of corrected reflectance for all 
clear-sky periods, starting with the first clear-sky composite period after 01-Aug. After the best-fit 
coefficients were determined, the new values were calculated using the polynomial coefficients to 
replace contaminated pixels in each channel prior to the first clear pixel or after the last such pixel. 

Step 6: NDVI processing 

Because of imperfections in the bidirectional corrections of channels 1 and 2, the NDVI values 
computed from atmospherically corrected NDVI were also retained. However, corrections for solar 
zenith angle were desirable hi view of the known dependence of the NDVI on the solar zenith angle. 
The coefficients of Sellers et al. (1994) were used for the various land cover classes. The new set of 
NDVI values was then computed for a reference solar zenith angle of 45 degrees based on the 
equations of Sellers et al. (1994). The NDVI values for the missing or contaminated pixels were 
interpolated as in Step 5 above. 

Step 7 : Compression 

BORIS staff copied the ASCII and compressed the binary files for release on CD-ROM. 


9.2.2 Processing Changes 
None. 

9.3 Calculations 

See Section 9.2.1. 

9.3.1 Special Corrections/Adjustments 

See Section 9.2.1. 

9.3.2 Calculated Variables 

See Section 9.2.1. 

9.4 Graphs and Plots 

None. 


10. Errors 


10.1 Sources of Error 

The sources of errors in LAI and FPAR are: (1) ground measurements, (2) use of algorithms based 
on Landsat TM that include image geometric and radiometric corrections and atmospheric correction, 
(3) error in determining the response difference between TM and AVHRR sensors, (4) the error 
associated with the AVHRR NDVI product. 

The major sources of error in AVHRR NDVI are due to the inaccurate knowledge of atmospheric 
conditions during image acquisition (and thus the use of nominal values for atmospheric corrections) 
and imperfect modeling of the bidirectional effects. 

10.2 Quality Assessment 

10.2.1 Data Validation by Source 

None given. 
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10.2.2 Confidence Level/Accuracy Judgment 

The FPAR images are more accurate than the LAI images. The relative error in LAI is estimated to 
be 20% for conifers and transitional forests; 25% for deciduous and mixed covers; and 30% for 
cropland, grassland, and tundra. The relative error in daily green FPAR is estimated to be less than 
10% for all cover types. 

10.2.3 Measurement Error for Parameters 

See Section 10.2.2. 

10.2.4 Additional Quality Assessments 

None given. 

10.2.5 Data Verification by Data Center 

BORIS staff displayed the LAI and FPAR images as a visual check that the images were what was 
expected. 


11. Notes 


11.1 Limitations of the Data 

The calculated LAI and FPAR maps should be considered the best estimate for this moderate 
resolution. The distribution patterns and the magnitude of LAI and FPAR values for each cover type 
are indeed very reasonable. The resampling scheme used for the original channel images makes the 
actual resolution coarser than 1 km. It is on the order of 2-3 km in reality (Cihlar et ah, 1997b). 
Subpixel water bodies and mixture of different cover types in a pixel can incur errors in the calculated 
results. Lack of data makes the calculations for rangeland, tundra, and mixed wood types least reliable. 

11.2 Known Problems with the Data 

None. 

11.3 Usage Guidance 

Before uncompressing the Gzip files on CD-ROM, be sure that you have enough disk space to 
hold the uncompressed data files. Then use the appropriate decompression program provided on the 
CD-ROM for your specific system. 

11.4 Other Relevant Information 

None. 


None given. 


12. Application of the Data Set 


13. Future Modifications and Plans 


None given. 
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14. Software 


14.1 Software Description 

A Fortran program, "lcc_1200.f", is provided for the calculation of longitude and latitude for any 
given pixel and line in the image. When using the software, you need to input the pixel and line 
numbers relative to the top-left pixel, which is (1,1). Gzip (GNU zip) uses the Lempel-Ziv algorithm 
(Welch, 1994) used in the zip and PKZIP commands. 

14.2 Software Access 

A version of the lcc_1200 code written in C is available on request from contacts listed in Section 
2.3. Gzip is available from many Web sites across the Internet (for example, FTP site 
prep. ai.mit.edu/pub/gnu/gzip-*. *) for a variety of operating systems in both executable and source 
code form. Versions of the decompression software for various systems are included on the 
CD-ROMs. 


15. Data Access 

The RSS-07 regional LAI and FPAR images are available from the Earth Observing System Data and 
hiformation System (EOSDIS) Oak Ridge National Laboratory (ORNL) Distributed Active Archive 
Center (DAAC). 

15.1 Contact Information 

For BOREAS data and documentation please contact: 

ORNL DAAC User Services 
Oak Ridge National Laboratory 
P.O. Box 2008 MS-6407 
Oak Ridge, TN 37831-6407 
Phone: (423) 241-3952 
Fax: (423) 574-4665 

E-mail: omldaac@oml.gov or ornl@eos.nasa.gov 

15.2 Data Center Identification 

Earth Observing System Data and Information System (EOSDIS) Oak Ridge National Laboratory 
(ORNL) Distributed Active Archive Center (DAAC) for Biogeochemical Dynamics 
http://www-eosdis.oml.gov/. 

15.3 Procedures for Obtaining Data 

Users may obtain data directly through the ORNL DAAC online search and order system 
[http://www-eosdis.oml.gov/] and the anonymous FTP site [ftp://www-eosdis.oml.gov/data/] or by 
contacting User Services by electronic mail, telephone, fax, letter, or personal visit using the contact 
information in Section 15.1. 

15.4 Data Center Status/Plans 

The ORNL DAAC is the primary source for BOREAS field measurement, image, GIS, and 
hardcopy data products. The BOREAS CD-ROM and data referenced or listed in inventories on the 
CD-ROM are available from the ORNL DAAC. 
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16. Output Products and Availability 


16.1 Tape Products 

The data can be made available on 8-rnrn or Digital Archive Tape (DAT) media. 

16.2 Film Products 

None. 

16.3 Other Products 

These data are available on the BOREAS CD-ROM series. 
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None. 


18. Glossary of Terms 


None given. 


ABC 3 

AEAC 

APT 

ASCII 

AVHRR 

BOREAS 

BORIS 

BP I 

BSQ 

CCRS 

CCT 

CD-ROM 

CECANT 

DAAC 

DAT 

DN 

EOS 

EOSDIS 

EROS 

FASIR 

FPAR 

GAC 

GEOCOMP 

GIS 

GSFC 

HRPT 

IFC 

IFOV 

ISLSCP 

LAC 

LAI 

LCC 

MRSC 

MSR 

NAD 8 3 

NASA 


19. List of Acronyms 

Atmosphere, Bidirectional and Contamination Corrections of CCRS 
Albers Equal-Area Conic 
Automatic Picture Transmission 

American Standard Code for Information Interchange 

Advanced Very High Resolution Radiometer 

BOReal Ecosystem-Atmosphere Study 

BOREAS Information System 

Bytes Per Inch 

Band Sequential 

Canada Centre for Remote Sensing 
Computer-Compatible Tape 
Compact Disk-Read-Only Memory 

Cloud Elimination from Composite Using Albedo and NDVI Trend 

Distributed Active Archive Center 

Digital Archive Tape 

Digital Number 

Earth Observing System 

EOS Data and Information System 

Earth Resources Observation System 

Fourier-Ad justment , Solar Zenith Angle corrected, Interpolated, 
Reconstructed 

Fraction of PAR absorbed by plant canopies 
Global Area Coverage 
Geocoding and Compositing System 
Geographic Information System 
Goddard Space Flight Center 
High-Resolution Picture Transmission 
Intensive Field Campaign 
Instantaneous Field-of -View 

International Satellite Land Surface Climatology Project 

Local Area Coverage 

Leaf Area Index 

Lambert Conformal Conic 

Manitoba Remote Sensing Centre 

Modified SR 

North American Datum of 1983 

National Aeronautics and Space Administration 
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NDVI 

NEdT 

NOAA 

NRL 

NSA 

ORNL 

PAR 

PASS 

PRT 

RSS 

SMAC 

SR 

SSA 

SST 

TIROS 

TM 

TOA 

TRAC 

URL 


- Normalized Difference Vegetation Index 

- Noise Equivalent Differential Temperature 

- National Oceanic and Atmospheric Administration 

- Naval Research Laboratory 

- Northern Study Area 

- Oak Ridge National Laboratory 

- Photosynthetically Active Radiation 

- Prince Albert Satellite Station 

- Platinum Resistor Thermometer 

- Remote Sensing Science 

- Simplified Method for Atmospheric Correction 

- Simple Ratio 

- Southern Study Area 

- Sea Surface Temperature 

- Television and Infrared Observation Satellite 

- Thematic Mapper 

- Top-of -the-Atmosphere 

- Tracing Radiation and Architecture of Canopies (a LAI and FPAR 
optical instrument) 

- Uniform Resource Locator 


20. Document Information 

20.1 Document Revision Date 

Written: 25-Sep-1997 
Last Updated: 26-Jul-1999 

20.2 Document Review Date(s) 

BORIS Review: 06-Oct-97 
Science Review: 10-Jan-98 

20.3 Document ID 


20.4 Citation 
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products by the MRSC in Winnipeg, Manitoba. The present level-4c product was created by CCRS 
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If using data from the BOREAS CD-ROM series, also reference the data as: 

Chen, I.M. and I. Cililar, "Retrieval of Boreal Forest Leaf Area Index From Multiple Scale 
Remotely Sensed Vegetation Indices." hi Collected Data of The Boreal Ecosystem-Atmosphere Study. 
Eds. I. Newcomer, D. Landis, S. Conrad, S. Curd, K. Huemmrich, D. Knapp, A. Morrell, I. 
Nickeson, A. Papagno, D. Rinker, R. Strub, T. Twine, F. Hall, and P. Sellers. CD-ROM. NASA, 
2000 . 

Also, cite the BOREAS CD-ROM set as: 

Newcomer, I., D. Landis, S. Conrad, S. Curd, K. Huemmrich, D. Knapp, A. Morrell, I. 
Nickeson, A. Papagno, D. Rinker, R. Strub, T. Twine, F. Hall, and P. Sellers, eds. Collected Data of 
The Boreal Ecosystem-Atmosphere Study. NASA. CD-ROM. NASA, 2000. 
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